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a^ 

j>^, Abstract 

^: 

We develop the idea that, in quantum gravity where the horizon fluctuates, 
a black hole should have a discrete mass spectrum with concomitant line 
emission. Simple arguments fix the spacing of the lines, which should be broad 
, but unblended. Assuming uniformity of the matrix elements for quantum 

CNj \ transitions between near levels, we work out the probabilities for the emission 

of a specified series of quanta and the intensities of the spectral lines. The 
• thermal character of the radiation is entirely due to the degeneracy of the 

, levels, the same degeneracy that becomes manifest as black hole entropy. 

Q>^ I One prediction is that there should be no lines with wavelength of order the 

I black hole size or larger. This makes it possible to test quantum gravity with 

O^' black holes well above Planck scale. 
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In Hawking's conception of the origin of black hole radiance, modes prepared in the 
pure vacuum state at X_ contract, pass through the center of the collapsing star, and make 
it out to X+ after spending a long time "skimming the horizon" . The mode frequenies are 
redshifted down and the modes become populated by real quanta, the partners of those 
other members of the created pairs that have fallen through the horizon. The outgoing 
quanta have a thermal spectrum and statistics that bear the imprint of passage through 
the potential barrier surrounding the black hole. The universality of this thermal character 
is viewed as due to the long sojourn of the emergent quanta close to the horizon. 

Quantum fluctuations of the horizon, which are expected within any reasonable theory 
of quantum gravity, must change this simple picture. The Hawking quanta will not be 
able to "hover" at a nearly fixed distance from the horizon since the latter's very location 
has to fluctuate. Thus one suspects a modification of the character of the radiation when 
quantum gravity effects are properly taken into account, even for black holes very massive 
with respect to the Planck scale. In this letter we consider what modification of the naive 
Hawking radiance spectrum and statistics might reasonably be expected from quantum 
gravity. 

Quantum systems of finite size more often than not display a discrete energy spectrum. 
Since the dynamics of a black hole responsible for its unique character refer to the finite 
region enclosed by the horizon, one may thus expect the mass spectrum of a black hole to 
display discreteness. Long ago the case was made [0,^ that the black hole's horizon area 
should be quantized in integers. In particular (we use units with G = c = 1) 

A = ahn; n = l,2... (1) 

where a is a pure number. This type of quantization law has since been considered by many 
workers In particular, it is in agreement with the idea that area should be quantized 

in canonical quantum gravity The rule (|T]) implies that the mass of a nonrotating 
neutral black hole has a discrete spectrum. In general an energy level labeled by n will be 
degenerate; let its multiplicity be denoted by g{n). Further, it is natural to identify the 
entropy Sbh = \^/^ + const, of the black hole in the level n with \n.g{n). The rule (|l]) in 
conjunction with the natural assumption that g{l) = 1 or 5'bh(^ = 1) = (nondegenerate 
ground state) then forces us too choose g{n) = e"*^""^-*/'*. But g{n) must be integral, so 
a = 4\nk where k = 2,3,4 .. .. 

One of us ^ has argued that the choice = 2 is preferred because it would make the 
spacing in entropy of consecutive energy levels exactly one bit, an attractive value from 
the information-theoretic point of view. One convincing piece of evidence in favor of this 
choice is that it leads to g{n) = 2"~^ whereas the number of ways in which the black hole 
at level n can be built starting from a situation with A = (no black hole) and going up 
the staircase of levels by various combinations of steps is also 2"^^. Thus for k = 2 the 
degeneracy actually quantifyies the number of ways the black hole in the particular level 
could have formed. Put another way, the decay from the n-th level to the no-black hole 
state in such scheme can be accomplished in exactly 2""^ sequences of steps [0. There is 
thus weighty evidence in favor of = 2 or a = 41n2. 

With this choice the energy spacing between consecutive levels for M ^ h corresponds 
to the (fundamental) frequency 
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Clearly the radiation emitted as the black hole decays will be concentrated in lines at in- 
teger multiples of the fundamental frequency w. Even if some broadening occurs or if the 
spectrum is blurred by two-quanta emission per jump, it is clear that the spectrum will be 
radically different from the one entertained in the standard discussion of Hawking radiance. 
One consequence of the physics here considered is that no radiation (or very little) should 
be emitted below w. Observation of this type of spectrum for any black hole would immedi- 
ately make quantum gravity effects observable well above the Planck scale. Herein lies the 
importance of our considerations. 

Information about the intensity of the hnes emitted by the black hole in the various 
particles available, and the statistics of these emissions can be obtained by first focusing on 
the de-excitation probabilities. We define these as follows. Describe the decay of the black 
hole during any interval of observer time At by giving a sequence of integers {ni, n2, ■ ■ ■ , n^}. 
This means that during At the black hole first jumped down ni elementary levels in one go, 
then 712 levels, • • •, and concluded by jumping down rij levels. In the process it emitted a 
quantum of some species of energy riihzu, then a quantum of energy n2hzu, and so on. 

These sequences can occur with any length j. The sequence of zero length, formally 
written {0}, represents the eventuality that the black hole did not decay at all in the 
interval At. One can associate with any sequence of length j a conditional probability 
PAt({^i, ^2, • • • , nj}\j) whose dependence on the species of quanta we ignore in this prelim- 
inary exploration. Obviously normalization requires that 

PAt{{ni,n2,---,nj}\j) = 1, (3) 

{ni;n2,---,nj} 

where the all the are nonvanishing. 

Even before dealing with the conditional probability PAt{{f^iTf^2, " " " given that 

the hole decayed by a sequence of j jumps, we must evaluate the probability PAt{j) for ^ 
sequence of exactly j jumps (and j emitted quanta) in time At. If At is not vastly large, 
the black hole will have changed very little after time At in terms of its emission properties 
(it will have nearly the same mass). Then the probabilities PAtU) for the next time interval 
At can be taken as equal to the previous PAt(j)- On the basis of this and the fact that one 
does not attach significance to a failure of the hole to jump at some point in a sequence, we 
can write, for instance, 

P2At(l) = PAt{0)pAt{l) + PAt(l)pAt(0) = 2pAt(0)pAt(l) (4) 

In words: a jump by one level in the interval 2At can take place either by the jump taking 
place in the first half of the interval with no jump in the second half, or by no jump taking 
place in the first half with the jump by one step taking place in the second half. The above 
equality may be generalized to any odd j by 

P2At{j) = 2pAt(0)pAt(j) + 2pAt{l)pAt{j -!) + ••• + 2pAt(j72 + l/2)pAt(j72 - 1/2) (5) 

and to any even j by 
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P2U3) = 2pAt{0)pAt{j) + 2pAt{l)pAt{j - 1) + ■ ■ ■ + [pAtUmY (6) 



These equations may be supplemented by the very obvious one 

P2At(0) = [pAmf (7) 

which just says that the survival probability of the black hole in a given level over time 2At 
is product of the survival probability over time At by itself. The solution of Eq. (|^ is 

PAtiO) = e-^*/- (8) 

where r is a survival timescale to be determined. Substituting Eq. (||) into Eq. gives a 
simple functional equation for pAt(l) whose general solution is 

PA,(l) = (At/r*)e-^*/- (9) 

where r* is a positive constant — possibly distinct from r — to be determined. Now sub- 
stituting Eqs. into Eq. (|) gives a functional equation for PAt(2) which we solve 
by 

PAi(2) = (l/2)(At/r*)2e-^*/- (10) 

On the basis of Eqs. (|])-(|r3) l^t us guess that for any j 

PAt(j) = (l/j!)(AVr*ye-^*/- (11) 

To prove this guess for any odd j we substitute Eq. (H) and our guess in Eq. (^), factor out 
the common factor (At/r*)-'e~-'^*/'^, and multiply the resulting equation by j\ to get 

9.-1 1 . , , j(j-l)(j-2) , , j(j-l)- •■072 + 1/2) 

But this is plainly correct by the binomial expansion of (1 + 1)-^. Thus Eq. ([TT| ) is correct 
for all odd j. The proof of Eq. (|n|) for even j proceeds similarly from Eq. (^). We must 
now check the normalization of pAt{j)' 

EPAi(j)=exp(At/r*-At/r) (13) 
j 

Obviously we must set r* = r to get normalization. 

Thus we have shown that the lengths of the sequence {rii, n2, ■ ■ ■ , n^} (in other words, the 
number of quanta emitted during At) follows a Poisson probability distribution, [Eq. (pUf) 
with r* = r] depending on a single parameter r. It is a consequence of Eq. ([TT|) that the 
mean value of j is simply (j) = At/r. 

Let us now return to the discussion preceding the introduction of pAt{j)- Given that the 
black hole has decayed during the interval At by a sequence {ni,n2,- ■■,nj} of length j, 
what is the corresponding conditional probability distribution PAt{{ni, n2, - ■ ■ , nj}\j) ? For 
j = 1 the answer is not difficult to guess. When the black hole jumps rii steps down the 
ladder of levels, the degeneracy g{n) changes by a factor of 2~"^ If we assume the "matrix 
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element" for a jump is the same for all ni 7^ 0, then the probability for the jump must be 
proportional to the final level's degeneracy. Thus 

PAi({ni}|l)(x2-"^ (14) 

(The case ni = is not included here because it belongs to j = 0). From the normalization 
J2'^^=i PAt{{ni}\l) = 1 we see immediately that the proportionality constant in Eq. (|T3[) is 
unity. Combining Eq. (|l^) with Eq. (|D we get for the a priori probability that the hole 
jumps once rii steps in At 

PAtiini}) = PAti{ni}\l)pAtil) = (At/2"V)e-^*/^ (15) 

Now consider the case j > 1. Up to now we have not regarded At/r as having any 
particular value. However, it is clear from Eq. (|Tl|) that when At ^ r, pAt{j) is very small 
for j > 1. We are thus interested in calculating PAt{{ni,n2, ■ ■ ■ ,nj}\j) only for At/r not 
small compared to unity. However, it is a useful strategy to imagine the interval At divided 
into a large number N ^ j of sub intervals of equal durations e = At/N. Because quantum 
transitions "do not take time", we can think of the j jumps as taking place each in one of 
j of the subintervals. There is no need to consider two or more jumps in one subinterval 
because, by our earlier remark, the probability of such an event would be very small because 

is large. 

Now the j "active" subintervals can be chosen out of the subintervals in a total of 
N\/j\{N — j)! ~ / j\ ways, where the approximation is a good one because j <^ A^. Thus 
the analog of Eq. ([T5|) is 

PMnu n2,---, n,}) ^ {N^ /j\)mn,}\l)P,{{n2}\l) ■ ■ ■ mn,}\l)[p,{l)y[p,{0)f-^ (16) 

where the last factor is the probability for the N — j subintervals during which the hole does 
not decay. Substituting from Eq. (lTl|) and replacing e ^ At/N we have 

PA*({ni, n2, ■ ■ ■ , n,}) = (1/j!) 2-("i+"^ +■■■+".) {At/rY e"^'!^ (17) 

Under the stated assumptions this equation should be exact as A^ ^ cxd. The set of prob- 
abilities is normalized in the sense that when summed for fixed j over all configurations 
{ni,n2,- ■ ■ -.nj) with every Uj starting from one, the result is exactly PAt{j), so that the 
further sum over j yields unity. Another thing to notice is that 

PAtiini, ^2, ■ ■ ■ , n,}\j) = PAt{{nun2, n,})/pAtU) = 2-("^+-^+-+-^). (18) 

Thus given that the black hole jumped j times in At, the distribution of the sizes of the 
jumps is independent of At as well as any scale set by the black hole. (In fact, we could obtain 
Eq. (|18D directly by repeating our argument that if the transition matrix element is constant, 
the probabilities depend only on the multiplicity of the final level.) One consequence of this 
is that the spectrum and the statistics of the radiation are steady so long as the mass of the 
hole changes by little. 

Let us now estimate the parameter r. The conditional probabilities Eq. ([I8|) are nor- 
malized. Since (rij) = J2'^i=i ni2~^^ = 2, each one of the j quanta carries energy 2hw on 



5 



average. Since the mean value of j is At/r, we see that the mean decrease in black hole 
mass during At is just 2hwAt/T so that 



d{M)/dt = -2hwAt/T (19) 

Due to the thermal character of the radiation, we may estimate d{M)/dt by using the 
Stefan-Boltzmann law for a surface of area 167r(M)^ at temperature {h/87i){M)~^: 

d{M) ^ 7^ 
dt 15307r(M)2 ^ ' 

Here 7 accounts for two fact: the radiation comes out only in lines and not over the entire 
black-body spectrum, and the emission is by a variety of particles. The first fact tends to 
force 7 under unity, while the second promotes a value well above unity. On the balance 
we feel 7 should be regarded as being unity to within an order of magnitude. Inserting the 
value of w and comparing these equations gives us 



r 



38407"^ (M) In 2 (21) 



It is plain that r, the mean time between quantum leaps, is large compared to the dynamical 
black hole timescale (M). 

We now calculate the probability distribution p/^t{k\nk'Uj) that in the course of time At, 
the black hole emits k quanta of frequency n^ctj. If the hole makes j jumps in At, the number 
of ways of selecting the k quanta out of the j is Cl = i\[k\{i — fc)!]~^, and the probability 
of each such selection is 2~("i+"2H h"j-fc)(^2~"'*)*^ with none of the ni,n2, ■ ■ ■ ,nj_fc agreeing 
with rik- We must sum the product over all allowed values of rii 7^ n^, n2 7^ n^, etc.. Since 
Y^'^=i 2~" = 1, the sum over each of the j — k factors 2~"' gives a factor 1 — 2""'-'. Thus 

PAt{k\n,r^,j) = — ^ fl - (tz;^)' (22) 



k\{j -k)\\ 2"*/ V2"fe 

We now multiply this conditional probability by the absolute probability distribution PAt{j) 
in Eq. (pH]), and sum over all j > k. Since J2'jLs U — s)\ = , we have again a Poisson 
distribution: 

PM{k\nkw) = (1/A;!) (x^J^e'^-^ x„ = (At/r) 2"". (23) 

Is this result not at variance with ones conception of black hole radiation as thermal ? 
For thermal radiation the number of quanta in a given mode is distributed exponentially; 
here we have a Poisson distribution. However, the discrepancy is merely a result of the 
different sort of question we ask here: how is the number of quanta emitted during a time 
interval in all modes of frequency Uk^zr distributed ? 

To see that the Poisson distribution which answers this question is consistent with a ther- 
mal distribution, consider a blackbody cavity at temperature T. Let us look at all quanta 
contained in a random subvolume. Each quantum drawn from that subvolume is subject to a 
Boltzmann distribution; accordingly, the probability that we draw one quantum of frequency 
uj is A exp{—huj /T) where A is a normalization constant. Succesive drawings are independent 
drawings from the same distribution. Thus the probability that we draw exactly k quanta 
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of frequency Uk from a series of j drawings is C^A'' exp{—khuJk/T) y4-'~*T[^rf exp(— ^cUj/T) 
summed over all possible values of the frequencies Ui distinct from Uk- By normalization 
^ 6xp(— ^w/T) = 1, so if the ujk is missing from the sum, the corresponding factor in 
the product, after the summation, will be 1 — A exp{—hu)k/T) . We get for k a probability 
distribution exactly like Eq. (p2D with the replacement 2^"'- Aexp{—uJk/T). The per- 
fect analogy dispels the impression that the distribution for the black hole radiance is not 
thermal. 

Another interesting question is what is the probability that ki quanta are emitted into 
one group of modes, all of frequency UkZU, and /c2 into an equal number of modes of the same 
frequency (distinguished from the first group, say, by angular momentum) ? Since there are 
C-kl^^^ ways to divide ki + ^2 quanta of frequency Ukzu into one class with ki and one with 
k2, and a total of 2''^'^'^'^ such bipartitions, 

PAt{k,, k2\nkw) = pAt{ki + k2\nuw) 2-^=1+^^) C,'i+'=^ (24) 

With help of Eq. ( p3D this takes the form 

PAt{ki,k2\nk'Co) = PAt/2{ki\nkVj)pAt/2{k2\nk'Co) (25) 

That is, the probabilities for ki and k2 are independent and of Poisson type, but their means 
are, naturally, reduced to half of that for ki + ^2. The generalization to the probability for 
many sets of modes of equal frequencies is obvious. 

Let us now inquire into the probability PAt{k, l\nkW, riivj) for the emission of k quanta of 
frequency and / quanta of different frequency rii in time At. There are j![/c!/!(j — Z — A;)!]~^ 
ways for the k and / quanta to be chosen — without regard to identity — from the j. The 
probability of each such choice is 2~^'^'^^^'^^"'^^^-^\2~^^Y (2"""')'. We sum the product over 
all allowed values of rii ^ nk and rii ^ rii, etc. By the logic already invoked, the sum over 
Hi of each of the j — k — I factors 2~"' is replaced by a factor 1 — 2""'^ — 2~"'. Instead of 
Eq. (^ we get 

pUk. l\nkw, n,r.,j) = ^'^ _ (l ' ^ ' - l) - l) (26) 

We now multiply by PAt{j) of Eq. (|1T]), and sum over all j > A; + / to get the result 

PAt{k, lluk-cu, riiw) = PAt{k\nkw) PAt{l\niw) (27) 

Thus the probabilities for emission into the two frequencies are independent and of Poisson 
type. The generalization of this result to n frequencies is obvious. 

We conclude by commenting on the shape of the spectrum. The mean value of the number 
of quanta emitted at a given frequency uj = nw is, according to Eq. (^3]), proportional to 
At and to 2~". Since emissions at diverse frequencies are independent, it follows that 
the intensity of the line at frequency nvj must be proportional to n2~^ or equivalently to 
^g-87rMa;_ accordaucc with our earlier assumption, this is independent of particle species. 
The spectrum is to be compared with the Planckian spectrum (phase space factor aside) 
^^l^^huj/T _ 1). The black hole spectrum is obviously similar to that of thermal radiation of 
like temperature in a finite box. 



7 



We found earlier that the rate of emission of quanta is r~^. Thus the width of the hues 
in frequency should be of order r~^. But we concluded that r is between a few hundreds 
and a few tens of thousands of M. We thus expect the width of the lines in frequency to be 
small compared to M~^, which is well below the size of the frequencies, Eq. Therefore, 
the various black hole lines are unlikely to overlap; the black hole spectrum must be clearly 
discrete. This conclusion differs from Hawking's original prediction even for massive black 
holes. Fluctuations of the horizon are here seen as manifested in a discretization of the 
radiance's spectrum. No quanta are expected with frequencies below . These predictions 
will become amenable to experimental check if primordial black holes are ever found. They 
open a window upon quantum gravity that does not require reaching down to Planck scale 
physics. 

V. M. is grateful to the Racah Institute, Hebrew University, for hospitality. J. D. B.'s 
work is partly supported by a grant from the Israel Science Foundation, administered by 
the Israel Academy of Sciences and Humanities. 



8 



REFERENCES 



[1] S. W. Hawking, Commun. Math. Phys. 43, 199 (1975). 

[2] L. Parker, Phys. Rev.D12, 1519 (1975); J. D. Bekenstein, Phys. Rev.D12, 3077 (1975); 

R. M. Wald, Commun. Math. Phys. 43, 199 (1975). 
[3] J. D. Bekenstein, Lett. Nuovo Cimento 11, 467 (1974). 

[4] V. Mukhanov, Pis. Eksp. Teor. Fiz. 44, 50 (1986) [JETP Letters 44, 63 (1986)], and 
in Complexity, Entropy and the Physics of Information, SFI Studies in the Sciences of 
Complexity, vol III, ed. W. H. Zurek (Addison- Wesley, New York 1990). 

[5] Yu. I. Kogan, Pis. Eksp. Teor. Fiz. 44, 209 (1986) [JETP Letters 44, 267 (1986)]; M. 
Schiffer, "Black Hole Spectroscopy" (unpublished); J. Bellido, "Quantum Black Holes", 
report |hep-th/9302T2^ (unpublished); M. Maggiore, Nucl. Phys. B429, 205 (1994); Y. 
Peleg, report [hep-th 9307057] (unpublished); C. Lousto, (unpublished). 

[6] U. H. Danielsson and M. Schiffer, Phys. Rev. D48, 4779 (1993). 

[7] A. Ashtekar, C. RoveUi and L. Smohn, Phys Rev Letters 69, 237 (1992). 



9 



